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We demonstrate quantum interference between indistinguishable photons emitted by two nitrogen- 
vacancy (NV) centers in distinct diamond samples separated by two meters. Macroscopic solid 
immersion lenses are used to enhance photon collection efficiency. Quantum interference is verified 
by measuring a value of the second-order cross-correlation function (;'^'(0) = 0.35 ± 0.04 < 0.5. 
In addition, optical transition frequencies of two separated NV centers are tuned into resonance 
with each other by applying external electric fields. Extension of the present approach to generate 
entanglement of remote solid-state qubits is discussed. 



The interference of two identical photons impinging on 
a beamsplitter leads to perfect coalescence where both 
photons leave through the same output port. This fun- 
damental effect, known as Hong-Ou-Mandel (HOM) in- 
terference [ij, is a consequence of bosonic statistics for 
indistinguishable particles. HOM interference has been 
demonstrated using single photon pairs from paramet- 
ric down-conversion |2| and delayed photons from a sin- 
gle photon source [3|-l5|- HOM interference has recently 
drawn attention as a resource for entanglement genera- 
tion between distinct single-photon emitters with many 
potential applications in quantum information science l(3| . 
The effect has been observed for photons emitted by pairs 
of atoms [tI and trapped ions Q, and has been used 
for entanglement generation of remote trapped ions Q. 
While isolated atoms and ions, which are nominally iden- 
tical, are a natural source of indistinguishable photons, 
extending these ideas to condensed matter systems can 
be challenging since two solid-state emitters are generally 
distinguishable because of their different local environ- 
ments. This Letter demonstrates quantum interference 
of two photons produced by nitrogen-vacancy (NV) im- 
purities in distinct diamond samples separated by two 
meters. Complementin g th e recent work involving other 
solid-state systems |1CH13| . the present solid-state real- 
ization is particularly significant, since electronic and nu- 
clear spins associated with NV centers can be used as a 
robust solid-state qubit memory, yielding potential scal- 
able architectures for quantum networks |l4l . [l5| . Specifi- 
cally, in combination with a recent demonstration of spin- 
photon entanglement 16|, the present work paves the 
way for entanglement generation between remote solid- 
state qubits. 

Unlike those associated with atoms in free space, the 
optical properties of NV centers embedded in a solid state 
vary substantially from emitter to emitter, especially in 
distinct samples. This inhomogeneity is due to variation 
in the local environments of NV centers and, in particu- 
lar, to variation in the local strain. Furthermore, coinci- 
dence experiments are limited by the collection efficiency 



for light emitted by the NV center. While a wide variety 
of approaches are currently being explored to enhance the 
collection efficiency 17h20||. we here utilize solid immer- 



sion lenses (SILs) fabricated from bulk diamond [2l| to 
facilitate the efficient collection of narrowband photons 
with identical properties from distant diamond samples. 
The SILs improve the collection efficiency by minimiz- 
ing total internal refiection at the air-diamond interface, 
which is significant because of the high refractive index 
{ud = 2.4) of the diamond host. We measure enhance- 
ment factors in the range of 6-10, depending on the po- 
sition of the NV center inside the SIL. Very recently, 
microfabricated SILs have been used to observe HOM 
interference from two NV centers separated by roughly 
20 /im on the same diamond chip [10|. 

In our experiment, we use two 1.0-mm diameter SILs 
that are fabricated from bulk electronic grade diamond 
and cut along the (100) crystal plane. The SILs are 
kept at a temperature of 8K in continuous flow helium 
cryostats that are separated by two meters, as shown 
in Fig. [TJ We characterized the spectral properties of 
several NV centers using photoluminescence excitation 
(PLE) spectroscopy. This technique involves collecting 
the red-shifted FSB emission while an external-cavity 
diode laser is scanned across the ZFL transitions (see 
the inset of Fig. [1]). The resulting PLE spectra reveal 
FWHM linewidths in the range of 50-250 MHz for in- 
dividual transitions of NV centers in both SILs. These 
linewidths, although broadened by charge fluctuations in 
the vicinity of the NV center due to the ionization of 
nearby charge traps [221] , are comparable to the narrowest 
linewidths observed in both synthetic and natural bulk 
diamond samples [H, . 

To obtain identical photons from two NV centers, the 
NV centers need to have transitions that are spectrally 
overlapping, and the emission from these individual tran- 
sitions for each NV needs to be isolated. By performing 
simultaneous PLE scans on NV centers in the two SILs 
with a single laser, we can directly measure the relative 
detuning of their optical transitions. In our experimen- 
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FIG. 1. Schematic of the apparatus. Two diamond SILs con- 
taining multiple NV centers are housed in continuous helium 
flow cryostats 2 m apart. Each SIL is addressed by a separate 
confocal microscopy setup, which includes a path for excita- 
tion at 532 nm, a collection path for the phonon sideband 
(PSB), and a collection path for the zero-phonon line (ZPL). 
The three optical paths are superimposed using dichroic mir- 
rors. The ZPL collection path passes through a half-wave 
plate (HWP) and an additional band-pass filter (633-647 nm) 
before being coupled into a polarization-maintaining single- 
mode 50:50 fiber beamsplitter. The two beamsplitter out- 
put arms are connected to a pair of avalanche photodiodes 
(APDs), completing the Hanbury Brown and Twiss detection 
setup. An excitation laser at 637 nm is connected in place 
of ZPL APD 2 to acquire the photoluminescence excitation 
spectra. Electrodes for electric field tuning are installed in 
one cryostat. The inset shows a simplified level structure, 
including resonant excitation and emission into the ZPL at 
637 nm, red-shifted emission into the ground state PSB, and 
non-resonant excitation into the excited state PSB at 532 nm. 



tal sequence, a 5 //s pulse of green light initializes the 
NV center into the electronic spin sublevel of the triplet 
ground state with mg = (|0)) [24]. Therefore, we only 
collect fluorescence from the NV center when the laser is 
resonant with transitions from the |0) state to the |Ex} 
or |Ey) states, as shown in inset of Fig. [1] [2^. We need 
to select a pair of NV centers such that one transition 
in the first NV center is resonant, or can be tuned into 
resonance, with one transition in the second. 

We next demonstrate control over the optical proper- 
ties of NV pairs to compensate for strain-induced spectral 
inhomogeneities. We make use of the DC Stark effect to 
actively minimize the detuning between the selected tran- 
sitions 2^, 2^. Electric fields perpendicular to the NV 



axis vary the splitting between |Ex) and |Ey) states and 
parallel fields shift both transitions together. This allows 
complete control over the optical transition frequencies 



27l |. In order to apply the desired electric field, we 
place one of the SILs on top of a silicon wafer deposited 
with four electrodes, comprised of 40 nm Au on a Cr ad- 
hesion layer. The gate geometry is shown in Fig. [2ja). 



We apply a bias voltage, Kpp, on one of the electrodes 
while keeping the other three grounded. In Fig. the 
0) o |Ex) transition of NVl (blue) is tuned across the 
|0) |Ex} transition of NV2 (red) by varying the applied 
voltage Vapp from -30 V to 50 V, which creates an elec- 
tric field of up to 0.5 MV/m at NVl. At Vapp = -2.9 V, 
shown in Fig. [2jc), the detuning between the two tran- 
sitions is reduced to 25 ± 2 MHz from an initial value of 
270 ± 15 MHz. We measure hnewidths of 85 ± 2 MHz for 
NVl, which has been tuned, and 217 ± 4 MHz for NV2, 
which has not been tuned. Similarly, we do not observe a 
systematic change of the linewidths with applied external 
fields in several other NV centers. 

For the HOM interference measurement, we excite the 
NV centers with green light and collect the ZPL emis- 
sion. Because the green excitation ionizes charge traps 
in the diamond lattice (26j and these charge dynamics can 
change the total electric field at the NV center, the time 
during which we can collect fluorescence at the tuned fre- 
quency is significantly reduced. For this reason, we select 
NV centers for the HOM measurement whose transitions, 
shown in Fig. [2][d), are inherently detuned by 93 ± 15 
MHz and have linewidths of 88 ± 3 MHz and 106 ± 4 
MHz, eliminating the need for electric field tuning. 

We want to isolate the emission from the selected tran- 
sitions so that we only collect spectrally overlapping ZPL 
photons for the HOM measurement. ZPL emission is 
separated from the PSB emission using a dichroic mirror 
and a spectral band-pass filter, as described in the cap- 
tion to Fig. [TJ The linear and orthogonal polarization 
selection rules of the |0) ■f-J' |Ex) and |0) ■<-!> |Ey) tran- 
sitions allow us to select the emission from one of these 
transitions in each NV center by inserting linear polar- 
izers into the ZPL collection arms [2^ and setting the 
Objective HWPs, shown in Fig. [H to the correct angles. 
Because the ZPL collection used to measure the HOM 
interference and the resonant excitation used to perform 
the PLE scans follow the same optical path, we can use 
the PLE scans to set the correct polarization angle for 
the ZPL collection. Therefore, we can selectively collect 
photons emitted from the desired transitions under non- 
resonant excitation with green light. 

To confirm that we are addressing one single-photon 

emitter in each SIL, we infer the normalized, second- 

(2) 

order autocorrelation function gp^^ir) in a standard 
Hanbury Brown and Twiss setup by splitting the PSB 
emission in a 50:50 beamsplitter. We expect 5pgB(0) = 
for an ideal single-photon source, and the single-photon 
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FIG. 2. Electric field tuning of optical transitions, (a) Four 
Cr/Au gates deposited on silicon. The central gap has a di- 
ameter of 150 ^m. In this experiment, only one of the gate 
voltages was swept while the others were kept grounded, (b) 
PLE spectra for different applied gate voltages. The gate volt- 
age Vapp is varied from -30 to 50 V for NVl, which creates 
an electric field of up to 0.5 MV/m at NVl. The [0) o- |Ex) 
transition of NVl (blue) is tuned across the |0) |Ex) tran- 
sition of NV2 (red). For different Vapp, PLE spectra are offset 
by 20 kCts/s for clarity, (c) Linewidth measurement under 
electric field tuning. On resonance, the measured FWHM 
hnewidths are 85 ± 2 MHz for NVl (blue) and 217 ± 4 MHz 
for NV2 (red). The detuning of the optical transitions in two 
samples is 25 ± 2 MHz. (d) Linewidth measurement for the 
NV centers used for the HOM measurement. The measured 
linewidths are 88 ±3 MHz (green) and 106 ±4 MHz (orange), 
and the detuning is 93± 15 MHz without electric field tuning. 



nature of the emission is confirmed in Figs. [31[a,b). Res- 
onant photons from each NV center is sent to an individ- 
ual input port of a polarization-maintaining fiber-based 
beamsplitter. We balance the emission intensity by ad- 
justing the green excitation intensity for each NV center 
independently to obtain 1100 counts per second (Cts/s) 
per emitter at each output port of the beamsplitter. An 
additional HWP in one setup is used to adjust the po- 
larization matching of the photons at the beamsplitter. 
The output ports of the beamsplitter are connected to 
single photon detectors with timing resolution below 100 
ps. The cross-correlation between these detectors is eval- 
uated using a Time-Correlated Single Photon Counting 
Module with a resolution of 64 ps. 

Ideally, g'-^^O) = for a pair of indistinguishable pho- 
tons, but the minimal observable g^^-* (0) value increases 
in the presence of experimental noise, as described be- 
low. When the photons are distinguishable and the light 



FIG. 3. (a,b) Single emitter second-order autocorrelation 
functions, 

9psb(^)' of the PSB emission, inferred for the two 
NV centers used for the HOM measurement. The FWHM of 
the central antibunching features are 7.5 ± 0.1 ns for (a) and 
9.5 ± 0.2 ns for (b). (c,d) Demonstration of HOM interfer- 
ence from remote NV centers in the (c) distinguishable case 
and (d) indistinguishable case. The dashed lines indicate the 
limit expected from independent distinguishable single pho- 
ton sources at r = 0. Solid lines are a fit to the data based 
on the model described in the text. For the distinguishable 
case (c) , we find that the FWHM of the central antibunching 
feature is 9.2 ±0.4 ns and g'f\Q) = 0.54 ±0.04. For the indis- 
tinguishable case (d), we find that the FWHM of the central 
antibunching feature is 5.6 ±0.3 ns and g'^^\o) ~ 0.35 ±0.04. 
Error bars are estimated based on shot noise. 



intensity in both arms is balanced, the correlation mea- 
surement will yield 5(2' (0) = 0.5 fl. Measuring the 
cross-correlation function in the distinguishable case is 
equivalent to measuring the autocorrelation function of 
one emitter while the non-interfering emission from the 
other acts as uncorrelated noise, raising g^'^\0) from 
to 0.5. Therefore, a measurement of g^^\0) < 0.5 indi- 
cates quantum interference between photons emitted by 
the two single photon sources. 

Figures[3jc,d) show (?'-^-*(t) for two different settings of 
the HWP angle. In Fig. El^c), the angle is selected such 
that the emissions from the two NV centers are distin- 

('2) 

guished by their polarization, yielding g\^ {0) = 0.54 ± 
0.04. In Fig. [IJd), the photons are indistinguishable 
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when their polarizations are parallel, yielding g^^^ '{0) = 
0.35 ± 0.04. In terms of the visibility of the HOM inter- 
ference, defined as r/ = [gf\o) — g^^\o)] / gf\o) , we find 
rj — 35 ± 9%. This ?7 > clearly demonstrates quantum 
interference between photons emitted by two NV centers 
separated by 2 m. 

We next turn to the detailed discussion of our exper- 
imental observations. We first consider the sources of 
noise that will cause our result to deviate from the ideal 

(2) 

case (0) = 0. The APD dark counts and fluores- 
cence background from our samples will lead to coinci- 
dence events, independent of the emission from the NV 
centers. Background light and the dark counts of our de- 
tectors contribute 80 Cts/s out of the total 1100 Cts/s 

(2) 

signal, raising g^^ ' (0) to 0.14. Because the NV center 
spin is not perfectly polarized under green illumination 
[2J] , we expect to collect emission from other transitions 
(e.g. 1^2) to \ms — ±1)); since this emission is assumed 
to be circularly polarized, it is only partially filtered by 
the polarizer. Because our collection objective has a large 
numerical aperture, the polarizations of the |0) |Ex} 
and |0) o |Ey) emissions are not perfectly orthogonal. 
Using PLE spectra acquired during the HOM data ac- 
quisition period, we measure that the selected transitions 
contribute a minimum of 94% of the total ZPL emission 

detected. Emission from other transitions at different fre- 

(2) 

quencies raises the value of g\^'iO) by 0.13. Finally the 
polarization-maintaining fiber-based beamsplitters intro- 
duce rotations to the polarization of the emission, which 

increases the distinguishability of the two photons. This 

(2) 

contribution raises the (0) value by 0.07. Consider- 
ing these factors, we expect experimental imperfections 

(2) 

to raise the (0) value to 0.34, which is in very good 
agreement with our experimental observations. 

We now analyze the temporal behavior of the inter- 
ference data presented in Figs. |31[c,d). The cross cor- 
relation function at the output ports of a beam splitter 
whose input ports are balanced and driven by single pho- 
ton sources can be written as [ll| 

9^'Hr)^\g[l\r) + ^Ji{r) 
H{l-U\\r)gg\T)cosiAu;r)), (1) 

(2) 

where g^^ (r) are the second-order autocorrelation func- 
tions for each single photon source inferred from Figs. 
|3lja,b) and Auj is the frequency difference between the 
emitters. The slowly- varying component of the first- 
order correlation function for single photon sources is 
given by gll\T) ~ exp(— 7|t|/2) where 7 ~ 1/12 ns^-'- 
is the inverse lifetime of the emitter. In our model, we 
assume that the emission from the two NV centers is ra- 
diatively broadened with bandwidth 7, and that the 
center frequencies of the emitted photons are random and 
different for subsequent emissions. We assume the distri- 



bution of the center frequencies is given by the Lorentzian 
profile that we fit to the PLE spectra shown in Fig. [UJd). 
Two features are visible in Fig. Eld) for |t| < 10 ns, 
which correspond to different terms in Eq. ([T]). Quan- 
tum interference is described by the cosine term, which 
gives rise to a narrow interference feature, as described 
below. This feature sits on top of a broader antibunching 
feature given by the second-order autocorrelation func- 
tions. The amplitude of the interference term is fit using 
a phenomenological parameter ^. We convolve the quan- 
tum interference term with the frequency distribution, 
which washes out the cos (Awr) oscillations and deter- 
mines the 1/e full- width of the interference feature. We 
keep this width, which is determined independently by 
the PLE spectra to be 3.1 ns, constant for the fits in 
Figs. [3Ic,d) and find excellent agreement. 

The behavior of the measured g^^ (t) for |t| > 3.1/2 

(2) 

ns is determined solely by §1^ (t) . Using the model de- 
scribed in [28], we extract parameters, which are listed 
in the caption for Fig. [31 from the two-emitter cross- 
correlation datasets [Figs. [3Uc,d)] that are in good agree- 
ment with those extracted from the single-emitter auto- 
correlation datasets [Figs. [3lja,b)]. The small difference 
between the parameters of Figs. [3^a,b) and Figs. [3l^c,d) 
can be explained by drifts in laser intensity and focal spot 
position during the longer integration time used for the 
datasets shown in Figs. [3ljc,d). 

In summary, we have demonstrated the generation of 
indistinguishable photons from two spatially separated 
NV centers. Combined with the recent demonstration 
of entanglement between the electronic s pin of an NV 
center and the polarization of a photon [16[, our work 
paves the way for optically mediated generation of en- 
tanglement between remote solid-state quantum regis- 
ters. The techniques demonstrated here have yielded im- 
proved collection efficiency, control of the NV centers' op- 
tical transition frequencies via electric field tuning, and 
the ability to operate two independent setups simulta- 
neously over three days of continuous data acquisition. 
Additionally, implementing a resonant excitation scheme 
similar to that used in [l6| will likely result in narrow 
optical linewidths (23| . Such a scheme will also minimize 
the ionization of local charge traps, which increases mea- 
surement time while using electric field tuning with CW 
excitation at 532 nm. The important figure of merit for 
an entanglement experiment is the time required to gen- 
erate an entangled pair with fidelity greater than 50%. 
Using our currently available values for collection effi- 
ciency (4 X 10~^), narrow linewidths (50 MHz), and as- 
suming a repetition rate of 10*, we estimate that one 
entangled spin pair can be created within roughly ten 
seconds. Improved photon collection techniques that are 
currently being developed [2^ have the potential to in- 
crease this generation rate dramatically. Even with the 
currently estimated rates, though, the exceptionally long 
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nuclear spin memory times of NV centers [3C| may allow 
one to use such systems for the realization of solid-state, 
multi-node quantum networks. 
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